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1 Introduction. 

In the present note we give some "easily-stated" recurrences of a special type that generate the 
binary digits for some "complicated" real numbers, such as the one in the title. The binary digits 
of any real number t = (di^d^ ■ ■ ■ with 1 < t < 2 can be calculated by the formula 

d n = [t2 n ~ l \ -2[t2 n - 2 \, n>l. 

We here show a somewhat unexpected, "fancy" way to obtain the digits of some special multiples 
of v2, where it is possible to hide this calculation. 

This note is structured as follows. In Section [2] we first recall what is known about the Graham- 
Pollak sequence and its variants, which serve as the motivating examples for the definition of the 
so-called Graham-Pollak pairs in Section We then give the general theorem (Theorem 13. 2\i , 
which in particular provides a new extension of the original result due to Graham and Pollak. 
Finally, we give Corollary 3.3 as one surprising example of the general phenomenon. 

2 The Graham— Pollak sequence. 

As usual, denote by \x\ the greatest integer less than or equal to x € R, and by {x} the fractional 
part of x. Define the sequence (u n ) n >i by the recurrence 



u\ = 1, u. 



71+1 



\/2 (u n + ij 



n>l. (1) 



This sequence, which is also known as the Graham-Pollak sequence, first appeared in a proceedings 
paper of F. K. Hwang and S. Lin [8] in the framework of Ford and Johnson's sorting algorithm [4]. 
For the reader interested in the background of the algorithm, an updated exposition can be found 
in the third volume of D. E. Knuth's The Art of Computer Programming Ch. 5.3.1, pp. 188]. 
The sequence ([I]) was first investigated from a purely mathematical point of view by R. L. Graham 
and H. O. Pollak [6]. They found the particularly intriguing fact that 

d n = U 2n +l ~ 2U2n-l (2) 

gives the nth binary digit of ^2 = (1.011010100 . . .) 2 . 

This fact puzzled several authors since then, and it has often been included as a fun exercise 
in articles and books mostly on combinatorial number theory. We mention, for instance, P. Erdds 
and R. L. Graham [SJ p. 96], R. Guy [3 Ex. 30], R. L. Graham, D. E. Knuth, and O. Patashnik [SJ 
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Ex. 3.46]. More recent references are J. -P. Allouche and J. Shallit [H Ex. 45, p. 116] and J. Borwein 
and D. Bailey p} p. 62-63]. N. J. A. Sloane's online encyclopedia of integer sequences [12] gives eight 
sequences which are connected to the Graham-Pollak sequence ([1]), namely, A091522, A091523, 
A091524, A091525, A100671, A100673, A001521, and A004539. 

Recently [1 3|. E] , the present author found vast extensions of the Graham-Pollak sequence to 
parametric families of recurrences, where the initial value u\ = 1 is replaced by u\ = m and the 
in the recurrence is accordingly changed. However, the sequence is still wrapped in considerable 
mystery. Indeed, if we do not alter the \/2 in the recurrence, but on the other hand, allow only the 
1/2 to vary (if n is odd), some quite strange things happen: we get the digits of various different 
multiples of v2, whose digits are seemingly unrelated. We point out that if we let the 1/2 vary for 
n even instead (cf. |14^ Theorem 3.3]), such effects cannot be observed. 

3 Main Result. 

In this note we are concerned with the following type of recurrences. 
Definition 3.1. Let e £ K and define the sequence {v n ) n >\ by 



1. 



V n +l 



{[V2(v n + e)\, if n is odd; 
[s/2 (v n + |)J , if n is even. 



We call (e, t) a Graham-Pollak pair if the sequence 

d n = v 2n +i - 2v 2n -i, n>l, 
represents the binary digits of t; that is, t = (di-c^cfe ■ ■ -)2- 

Note that (1/2,^/2) is a Graham-Pollak pair according to the original result about the se- 
quence ([1]). Our main result is as follows: 

Theorem 3.2. A list of Graham-Pollak pairs is given by 



{(Ei,U): l<i<8}, 



where 



V2 i- 
1- -y- <ei < V2-1, U 

r- 19 r 

V2- 1 < e 2 < yV2- 13, t 2 
19 77 

yV2-13 <e 3 < yV2-54, t 3 

77 r- 309 r- 

— V2- 54 < £ 4 < — V^- 218, U 

^-218 < £ 5 < 12961 2 21 ° 37 72 - 916495974, t 5 = V~2, 

1296121037 ^ _ 916495974 ^ ^ < 79109 ^ _ ^ fa 759250125 

^V2 - 55938 <er<|>/2-3, tr 
^V2-3 < £8 < 



\/2- 1, 

8 8' 

32 32' 

128 128' 



72- 



314491699 



536870912 
46341 y/2- 19195 



536870912 



32768 

3 V2-i. 
2 2 



32768 ' 
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We first comment on a few aspects of the theorem. 



(a) A surprising feature of Theorem 13.21 is that as e varies continuously, the output makes dis- 
crete jumps among multiples of \[2. Figure 1 illustrates the various intervals for e and the 
corresponding numbers t appearing in Theorem | 



(b) The binary digits of y/2 are obtained for any choice of e in the interval 

[0.4959953 •• • , 0.5012400 ••• ). 
This slightly generalizes the original result of Graham and Pollak with e = 1/2. 

(c) There may well exist Graham-Pollak pairs besides those given in Theorem [321 However, it is 
easily checked with a computer that the range for admissible values of e cannot be too large. 
For example, for e = 0.2928 we get 6^067 = — 1, and for e = 0.7073 we have c?2293 = 2. Similar 
phenomena hold outside these bounds where the values d n = — 1 and d n = 2 are already 
obtained for smaller indices n. Therefore, possible new pairs can only arise in a very small 
neighborhood of e = 1 — y or e = y. 

(d) There is also a connection to normal numbers, which shows that a characterization result for 
Graham-Pollak pairs is very difficult to obtain. Suppose that there exists c\ G R such that 

{V2 < a < 1, for all k > 1. (3) 

Then - according to © below - the interval given for e\ can be enlarged to c\ ^1 — < 
£\ < v2 — 1. Similarly, if there is c 2 G M such that 

{3^ 2 k ~ 2 } > c 2 > 0, for all k > 1, (4) 

then the interval for eg can be enlarged to §\/2— 3 < £s < ^+c 2 ^1 — ^p^J ■ The inequality ((3|) 
implies that \[2 is not normal in base two, and ((4|) implies that 3\/2 is not normal [101 Ch. 1.8]. 

We conclude with a surprising example, which follows from Theorem 13.21 by the (rather plain) 
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observation that the number 1 — \ = 0.5086213 . . . lies in the interval given for eq. 



Corollary 3.3. Define the sequence (u> n )n>i by 
wi = 1, w n+ i 

Then for n > 31, W2n+i — 2u>2n~i is the (n + l)th binary digit of 759250125\/2- 



[V2 (w n + l- f^)J , ifnis odd; 
[y/2 (w n + ^)J , if n is even. 



4 Proof of Theorem 13.2 



First, let i G / := {1,2,... , 8} \ {5} and consider the pairs (e^tj) in the statement of Theorem 
Put 
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Figure 1: The sets {(e^, tj)} for i = 1, . . . , 8. 



with ai, f3i,k S Z and 2) = 1. It is easy to verify that a« + = 2 li+1 for i G /. Furthermore, let 
£i i and £2,1 be the endpoints of the associated interval for £j. We shall prove that for £1^ < £j < £2,1 
and /c > /j + 2 we have 

«2k = L^2 fc - 2 J +7i2 fc - /l ~ 2 , (5) 
^ + i = L^2 fc - 1 J+2 fc , (6) 

where 7^ = 2oti + This then implies that for k > U + 3, 

^fe+i - 2^-1 = L^2 fc - 1 j - 2[ta k - 2 \, 

which is the A;th binary digit of ij. In the final step we then show that formula © indeed holds 
true for < k < li + 1, which completes the proof. 

We first use induction to prove that if © holds for k = k + 2, then © and ([6]) hold for 
k> k + 2. Assume the validity of ©. We have to show that [y/2 (v 2 k + |)J = L^" 1 ] + 2 fc , which 
is equivalent to 

L^ 1 ] + 2 k < V2 (W~ 2 J + li 2 k ~ k - 2 + < L^2 fc - X J + 2 fe + 1, 
or in other words, 

< 2 k - l >- 1 ^ - + - 2 1 ^ + V2 La,^2 2 fc -^ 2 j - \a t ^2 2 fe -^ 1 j + ^ < 1. 

Since — Pi = 2a« and a« + A = 2 li+1 this is the same as 

< {ai V2 2 k - 1 '- 1 } - V2 { ai V2 2 k ~ l >~ 2 } + ^ < 1. (7) 
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Relation © is true since < {x} - V% {x/2} + y/2/2 < 1 for all 

Now, assume relation ©. We have to ensure that [y/2(v2k+i +e)J = [k2 k ~ 1 \ + j i 2 k ^ li ^ 1 , or 
equivalently, 

[t i 2 k ~ 1 \ +l i 2 k - h - 1 < V2(L^2 fc - 1 J +2 k + e) < L* i 2 fc - 1 J + 7 ;2 fc ^- 1 + 1. 
Here we end up with 

< {1 - V2){a,V2 2 k ~ h ' 1 } + V2 e < I, (8) 

which is true provided 1 — v2/2 < e < \[2j2. This interval includes all of the intervals £2,4) in 
Theorem 13.21 and hence there is no additional restriction on e. 

It remains to check the initial conditions. This task encompasses some straightforward calcula- 
tions; we only give the main steps. First, we have to guarantee that © is true for k = + 2. Of 
course, this crucially depends on the choice of e. Since v n (e) is non-decreasing for increasing values 
of e, there is at most one semi-open real interval [£i,i,£2,i) for e such that 

«2(i 4 +2) = L^2' l J + 7i = [ ai V2 - Pi\ + 2ai + # = [aiV2] + 2a;. (9) 

We will show that £2,4) = [Ci,i)C2,i)- It is not difficult to crank out a reasonable guess for £1^ 
with the help of a computer. In fact, ^2^+2) is a piecewise constant function in e with only a finite 
number of jump discontinuities. Thus, we can get a close approximation of £1^ by interval halving. 
Furthermore, from Definition 13.11 we see that £1^ (if it exists) has the form 5f v2 — dj for some 
integers Cj,dj E Z. We use Maple 11 (PolynomialTools [MinimalPolynomial] ) to calculate an 
approximate minimal polynomial of degree two with "small" coefficients to identify a conjectured 
value for £1^. Again, we have to ensure that the value still satisfies ([9|). 

As an illustration, let z = 6 and consider 

v 2 (h+2) = v e2 = v 62 {e), [aiV2\ + 2 ai = 2749487923. 

Figure 2 shows the location of the jumps in the graph of W62(£) for e G [0.40,0.60]. By the above 
procedure we find that £1^ is "close" to 

£1,6 = 1296121037^2/2 - 916495974 = 0.5012400 .... 

Once more, we use Maple with the ansatz e = — 5, where 5 denotes a small positive quantity, 
to show that e = £x,i is indeed the smallest value which satisfies ([9]). This is a symbolic computation 
and does not involve high-precision arithmetic. In a similar fashion, we show that £2,4 = It 
is important to note that the values of ^1,^3, . . . , i^^+ij+i remain unchanged for e S [£1,1,^2,1) for 
every fixed i £ I. Moreover, a routine calculation confirms that ([6]) is true for < k < li + 1. 

Finally, we have to treat the case i = 5, which is less involved than the cases i £ I. Here we 
directly show that 

V2k = [U2 k - 2 \ + 2 k ~ 2 and v 2k+1 = [t i 2 fc - 1 J + 2 k 

for k > 1, so that we do not have to bother about initial conditions. (We leave the details to the 
interested reader.) 

Summing up, we have that the intervals [£i t i, £1,2) are disjoint for i = 1, 2, . . . , 8 and completely 
cover [1 — ^). This finishes the proof of Theorem 13.21 □ 
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2,749,487,923- 

jumps of v[62] 



0.4 0.5 0.51 0.55 0.6 

e 

Figure 2: The jumps of v 62(e) for i = 6 and 0.4 < e < 0.6. 
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